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A technique called analytic perturbation theory, which respects the required analytic properties, consistent with causality, is applied 
to the definition of the running coupling in the timelike region, to the description of inclusive r-decay, to deep-inelastic scattering 
sum rules, and to the investigation of the renormalization scheme ambiguity. It is shown that in the region of a few GeV the 
results are rather different from those obtained in the ordinary perturbative description and are practically renormalization scheme 
independent. 


1 Analytic Running Coupling Constant 

The conventional renormalization-group resummation of 
perturbative series leads to unphysical singularities in the 
running coupling constant. For example, the usual QCD 
one-loop running coupling is f] 


The second, nonperturbative term, cancels the ghost 
pole. 

The above defines the running coupling in the space¬ 
like region. We, can also define a timelike (or s-channel) 
coupling a s (s)p which is related to the spacelike coupling 
through the following reciprocal relations: 


a FT (Q 2 ) 1 1 m 

4t r fa ln(Q 2 /A 2 ) ’ 1 1 

where /3o = 11 — 2n/3, n being the number of quark fla¬ 
vors. This evidently has a singularity (Landau pole) at 
Q 2 = A 2 , which is unphysical and inconsistent with the 
causality principle. Instead, we propose replacing per¬ 
turbation, theory (PT) by analytic perturbation theory 
(APT)0’0 to enforce the correct analytic properties, for 
example, that the running coupling be regular except for 
a branch cut for — Q 2 > 0, by adopting the dispersion 
relation 
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where the imaginary part is given by the perturbative 
result, that is 


p(a) = —Ima PT (—a — ie). (3) 

47T 

This leads to a consistent spacelike coupling, which in 
one-loop is: 
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where the contour in the first integral does not cross the 
cut on the positive real axis. In terms of the spectral 
function, then, 
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The spectral functions in 1-, 2-, and 3-loops are 
shown in Fig. [j], for n = 3. The areas under all these 
curves turn out to be the same, which implies a univer¬ 
sal infrared fixed point a *, 
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which is exact to all orders. at! It is also possible to prove 
that symmetrical behavior of the timelike and spacelike 
couplings is inconsistent with the required analytic prop¬ 
erties, that is, in any renormalization scheme 
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a We use the notation Q 2 = —q 2 , where Q 2 > 0 for spacelike mo¬ 
mentum transfer. 


a s (s)^a(Q 2 ), s = Q 2 . (9) 

This difference, which is important when the value of the 
running coupling is extracted from various experimental 
data,cl is demonstrated in Fig. || for the MS scheme. 
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Figure 1: The spectral functions: (a) and (b) are the one- and 
two-loop results; (c) is the two-loop approximation obtaineahas the 
first iteration of the exact renormalization-group equation; 121 (d) is 
the three-loop result in the MS scheme. 


Figure 2: The ratio of spacelike and timelike values of the running 
coupling constant in the APT approach, at one, two, and three 
loops. 


Because of the infrared fixed point, the running cou- 
pling implied by APT rises less rapidly for small Q 2 than 
does the usual perturbative running coupling. This is il¬ 
lustrated!] in Fig. [|. 

Another interesting fact is that within APT 
Schwinger’s conjecture □ about theuconnection between 
the P and spectral functions is valida Indeed, the /3 func¬ 
tion for the timelike coupling is 


_ s da s 
s 47t ds 



( 10 ) 


The above discussion assumes that the number of 
active flavors n is realistically small. However, if n is 
large enough, even the perturbative coupling can exhibit 
an infrared fixed point, at least in the two-loop level. 
This same fixed point defined by Eq. (|j) occurs in the 
analytic approach and we have 


— = I CA), n ~ 8 ’ 

4-7T \~Po/Pi, 9 < n < 16 . 


( 11 ) 


The transition between the nonperturbative and per¬ 
turbative fixed point occurs for n = 8.05, as is shown 
in Fig. |J. This is qualitatively consistent with the 
phase transition seen, for example, in nonperturbative 
approaches and lattice simulations^ 


2 Inclusive r Decay within APT 

The inclusive semileptonic r decay ratio for massless 
quarks is givenEl in terms of the electroweak factor Sew; 
the CKM matrix elements V u d , V us , and the QCD cor¬ 
rection A r , 

Rr = 35 E w(|Kd| 2 + |Ks| 2 )(1 + A t ) . (12) 



Figure 3: QCD evolution of the running coupling constants (de¬ 
fined in the spacelike region) compared to experimental dataQ Here 
a matching procedure across the various quark thresholds has been 
applied in the timelike regionjj 



Figure 4: The analytic infrared fixed point a* vs. n. 
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The correction A r may be written in terms of the func¬ 
tions r(s) or d(g 2 ), which are the QCD corrections to the 
imaginary part of the hadronic correlator II: Imll ~ 1+r 
and to the Adler D- function: D = — q 2 dH/dq 2 ~ 1 + d. 
The analytic properties of the Adler function allow us to 
write down the relations 

r°° d s 

d {<?) = -<? / 7 - W r(s), (13) 

i r s+ie dz 

r{s) = — / -d(z). (14) 

Because of the proper analytic properties (which are vio¬ 
lated in the usual perturbative approach) the QCD con¬ 
tribution to the t ratio is given by the two equivalent 
forms 



If one likes, one can think of r and d as effective running 
couplings in the timelike and spacelike regions, so as with 
the running couplings they may be expressed in terms of 
an effective spectral density, 
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possessing the same universal infrared limit as the run¬ 
ning coupling. 

In the ordinary perturbative approach, the function 
d(q 2 ) may be expanded in terms of the running coupling 
and in the third order is 

d PT (q 2 ) = a PT (-g 2 ) + d 1 [a PT (—g 2 )] 2 + d 2 [a PT (-g 2 )] 3 , 

(18) 

where we have introduced a = cxs/tt, and, numerically, 
for three active flavors, the coefficients are d^ s = 1.6398, 
d^ = 6.3710. Such is not the case in APT; rather, 
d ApT j s constructed from Eq. 0 ) with a spectral density 
obtained as the imaginary part of d PT on the physical 
cut: 

P eff (cr) = po(cr) +dipi(a) + d 2 p 2 (a), (19) 

where p„(cr) = Ima n+1 (—a — ie). We use the world aver¬ 
ages valuef] Ii T = 3.633 ±0.031. Our results are shown in 
Table [j], where for the sake of illustration we also show 
the PT results obtained by using the contour integral 
representation given in the second line of Eq. (ED- 

b This is consistent with the 1998 PDG valued which we extract 
a s R r = 3.650 ± 0.016. 


Table 1: The APT and PT parameters in the MS scheme extracted 
from r-decay. The first two rows refer to NNLO calculations, the 
last two rows to NLO calculations. The errors in the last digits are 
shown in parentheses. 


Method 

^ms (MeV) 

a(M 2 ) 

d(M 2 ) 

APT 

871(155) 

0.3962(298) 

0.1446(88) 

PT 

385(27) 

0.3371(141) 

0.1339(69) 

APT 

918(151) 

0.3983(236) 

0.1431(84) 

PT 

458(31) 

0.3544(157) 

0.1400(67) 


Most remarkably, APT exhibits very little renormal¬ 
ization scheme (RS) dependence. The r decay coeffi¬ 
cients di and d 2 are RS dependent, as are all but the 
first two beta function coefficients, defined by the renor¬ 
malization group equation 

= ~^-a 2 (l + aa +c 2 a 2 ). (20) 

d/A 2 

Because of the existence of RS invariants, one can inves¬ 
tigate the sensitivity of the predicted value of the QCD 
correction to the choice of RS by varying d\ and c 2 in 
a region where the degree of cancellation in the second 
RS invariant u> 2 = c 2 + d 2 — Cidi — d 2 does not exceed a 
specified limit, taking, for example, 


c 2 \ ± |d 2 1 + ci|di| ± d 2 ^ n 


( 21 ) 


That sensitivity is shown Edl in Fig. I based on = 
0.1881. Observe that the relative difference between the 
prediction of the lower corners of the domain defined by 
Eq. ( pl| ) is 0.8%, while in PT that difference is 5%. Note 
that the MS scheme lies outside this domain, as does 
the so-called V scheme, the prediction of which differs by 
only 0.2% from the MS value in APT, but by 66% in PT! 
To all intents and purposes, APT exhibits practically no 
RS dependence. (Because of the perturbative stability, 
the known three-loop level is quite adequate for these RS 
recalculations.) 


3 Deep-Inelastic Scattering Sum Rules 

At present, the polarized Bjorken and Gross-Llewellyn 
Smith deen inelastic scattering sum rules allow the 
possibility^^ as with r decay, of extracting the value of as 
from experimental data at low Q , here down to 1 GeV. 
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Figure 5: Contour plot of the APT correction A r at the three- 
loop order as a function of RS parameters d\ and C 2 . The dashed 
line indicates the boundary of the domain defined by Eq. 



Figure 6: F( n with 1-, 2-, and 3-loop QCD corrections vs. the 
coupling constant. 


3.1 Bjorken Sum Rule 

The Bjorken sum rule refers to the value of the integral of 
the difference between the polarized structure functions 
of the neutron and proton, 


/»00 

rr n =/ dx [g{{x,Q 2 ) - g^{x,Q 2 )] (22) 
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where the prefactor in the second line is the parton-level 
description. In the conventional approach, with massless 
quarks, the QCD correction is given by a power series 
similar to Eq. (|l^), with coefficients for three active fla¬ 
vors being = 3.5833, d^ s = 20.2153. However, this 
description violates required analytic properties of the 
structure function moments, which, as has been arguedp 
follow from the existence of the Deser-Gilbert-Sudarshan 
integral representation. 

Thus we adopt the analytic approach, which says 
instead 

A apt (Q 2 ) =6^\Q 2 ) + d 1 S^{Q 2 ) + d 2 S^(Q 2 ), (24) 
where 


This is not the case in the conventional PT approach, as 
is shownEd in Fig. |^. Second, the renormalization scheme 
dependence is again very small, so that various schemes 
which have the same degree of cancellation as the MS 
scheme give the same predictions all the way down to 
Q 2 ~ 1 GeV 2 . In contrast, it is impossible to make any 
reliable prediction for conventional PT, eazen if improved 
by the Pade approximant (PA) mnthodld for Q 2 below 
several GeV 2 . This is illustrated Ed in Fig. |l]. One can 
see that instead of RS unstable and rapidly changing PT 
functions, the APT predictions are slowly varying func¬ 
tions, which are practically RS independent. 


3.2 Gross-Llewellyn Smith Sum Rule 

A precisely similar analysis can be performed on the GLS 
sum rule, which refers to the integral 


'S'gls = 7 j 


dx 


Ff P (x, Q 2 ) + F^ p (x, Q 2 ) 


= 3 [l - Agls(Q 2 )] ■ 


(26) 

(27) 


Again, the APT approach leads to perturbative stability, 
and to practically no mnormalization scheme dependence 
down to very low Q 2 Ed 


1 f‘°° rlrr 

= JRT l [«St(-o - *)] • (25) 

that is, A APT is not a power series in a APT . 

Besides possessing the correct analyticity, the APT 
approach has two key properties in its favor. First, suc¬ 
cessive perturbative corrections are small, so that the 
two- and three-loop QCD results are nearly the same. 


4 Conclusions 

Our conclusions are four-fold. 

• APT maintains correct analytic properties (causal¬ 
ity), and allows for a consistent extrapolation be¬ 
tween timelike (r decay) and spacelike (DIS sum 
rules) data. 
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Figure 7: Renormalization scheme dependence of predictions for 
Abj vs. Q 2 for the APT and PT expansions. The solid curves, 
which are very close to each other, correspond to the APT result 
in the MS, A, B, and ECH schemes. The PT evolution in the MS, 
A, and B schemes are shown by dashed curves, the ECH scheme is 
indicated by a dotted curve, and the PA results in the MS, A, and 
B schemes are denoted by dash-dotted cimres. (The definitinm of 
the various schemes is given in our paperEJ) The SMC data|l3 is 
indicated Wi a square, the triangle is the E154 Antal 17 ! circles are 
E143 dataliij and the stars are recent E143 dataH 


• Three loop corrections are much smaller than in 
PT; thus, there is perturbative stability. 

• Renormalization scheme dependence is drastically 
reduced. The three-loop APT level is practically 
RS independent. 

• The values of A are larger in the APT approach 
than in the PT approach. Yet these values are con¬ 
sistent between timelike and spacelike processes, 
and consistent with the data. 

The work reported here is the beginning of a system¬ 
atic attempt to improve upon the results of perturbation 
theory in QCD. In the future we will treat in detail the 
significance of power corrections which come from the 
operator product expansion, and examine the effect of fi¬ 
nite mass corrections, which necessitate a more elaborate 
analytic structure. 
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